Transverse power flow reversing of confined waves in extreme nonlinear metamaterials 
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We theoretically prove that electromagnetic beams propagating through a nonlinear cubic meta- 
material can exhibit a power flow whose direction reverses its sign along the transverse profile. This 
effect is peculiar of the hitherto unexplored extreme nonlinear regime where the nonlinear response 
is comparable or even greater than the linear contribution, a condition achievable even at relatively 
small intensities. We propose a possible metamaterial structure able to support the extreme con- 
ditions where the polarization cubic nonlinear contribution does not act as a mere perturbation of 
the linear part. 

PACS numbers: 78.67.Pt, 42.65.Tg 



(N 



CZ3 . 

O . 
> ' 

or 
q : 

o ■ 

<Z3 ■ 

>v 

^ : 



> 
O 



Os 
O 



X 



The ability of manufacturing metamaterials with pre- 
scribed and anomalous values of permittivity e and per- 
meability /i has triggered an intense research effort aimed 
at investigating novel regimes of linear electromagnetic 
propagation and suitable configurations have been de- 
vised for observing remarkable effects such as, for exam- 
ple, superlensing |l|, |2j, optical cloaking [3|, |j|, guidingof 
nanometric optical beams [j| and photonic circuits • 
In the nonlinear realm, the nonlinear properties of left- 
handed metamaterials have been investigated [8( together 
with various soliton manifestations [t| 03, [ll[ . Propaga- 
tion in metamaterials exhibi ting cubic nonlinear response 
has also been considered [T^.ll3] and, for ultra-short pulse 
nonlinear dynamics, it has been suggested that metama- 
terial linear property tailoring allows the observation of 
different nonlinear regimes [14| . 

In this Letter we show that a metamaterial with a very 
small linear dielectric constant and exhibiting a nonlin- 
ear cubic response is able to support nonlinear guided 
waves whose Poynting vector has the very peculiar prop- 
erty of being parallel and anti-parallel to the propaga- 
tion direction in different transverse portion of the field. 
This novel phenomenology is a consequence of the fact 
that, since the metamaterial linear dielectric permittiv- 
ity can be arbitrary small, the nonlinear contribution to 
the dielectric response can easily (i.e. at low intensities) 
be made comparable or greater than the linear part so 
that, the sign of the overall dielectric response can be 
different for different intensities. In the presence of an 
electromagnetic beam this implies that conditions can be 
found so that the effective dielectric response has differ- 
ent signs on the propagation axis and at its lateral sides. 
Therefore the transverse reversing of the power flow is 
understood since, for a monochromatic Transverse Mag- 
netic (TM) field mainly propagating along a given di- 
rection, the Poynting vector globally lies along the same 
mean propagation direction and its sign coincides with 
that of the total effective dielectric constant. In order 
to discuss this effect on a feasible situation, we consider 
TM electromagnetic propagation in a defocusing nonlin- 



ear cubic metamaterial and we analytically obtain a class 
of nonlinear guided waves exhibiting the aforementioned 
transverse power flow reversing. It is remarkable that 
the power flow reversing effect can be observed even at 
very low intensities since it is a consequence of the in- 
terplay between the linear and nonlinear contributions 
to the dielectric response, the former being very small in 
the considered metamaterials and the latter being pro- 
portional to the intensity. To the best of our knowl- 
edge, the transverse power flow reversing predicted in 
the present Letter is the first example of effects charac- 
terizing a novel extreme nonlinear regime where the non- 
linear contribution to the medium polarization does not 
play the role of a mere perturbation to the linear part. 
The question naturally arises as to whether a medium 
exists or can be conceived where the range of electro- 
magnetic intensities, for which its nonlinear response is 
purely cubic, is so large to produce a huge nonlinear re- 
sponse. At first one may reject this possibility since the 
cubic nonlinear response generally arises from a pertur- 
bative description of radiation-matter interaction so that 
the nonlinear polarization necessarily is a small correc- 
tion to the linear part. However, exploiting the avail- 
ability of metamaterials with somehow prescribed val- 
ues of the dielectric permittivity, we propose that in a 
suitable sub-wavelength layered structure, consisting of 
alternating slabs of a metamaterial with negative dielec- 
tric constant and a standard nonlinear cubic medium, 
the effective electromagnetic response is purely cubic in 
an intensity range where the nonlinear cubic term can 
exceed the linear contribution. 

Consider a monochromatic electromagnetic field 
(whose time variation is assumed to be e~ llAjt , where u> is 
the angular frequency) propagating through a nonlinear 
metamaterial characterized by the constitutive relations 
(holding for the field complex amplitudes) 

D = e eE-e x[(E-E*)E + 7(E-E)E*], 

B = [iofxH, (1) 

where e > and \x > are the linear permittivity and 
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FIG. 1: Nonlinear guided waves transverse profile of u x (panel 
(a)) and of u z (panel (b)) at different values of u xoo in the 
range of Eq.©, for 7 = 0.5. 



permeability, respectively, whereas \ > and < 7 < 1 
are the parameters characterizing the cubic defocusing 
nonlinear response. We focus our attention on trans- 
verse magnetic (TM) nonlinear guided waves propagating 
along the z— axis of the form 



E(x,z) 



H(x, z) = e 




/MO 



e y (2) 



where £ = ^/ejl(uj/c)x, ( = ^/ejl(uj/c)z (c is the speed of 
light in vacuum) are dimensionless spatial coordinates, 
P is a real dimensionless propagation constant and u x 
and u z are dimensionless electric field components. Sub- 
stituting the fields of Eqs.© into Maxwell equations 
VxE = iijjH and V x H = — iujT) and using the consti- 
tutive relations of Eqs.|T|) we get 

= [(/3 2 -l) + (l+7K + (l-7)^K, 
^ - ^ = [-1 + (1 - 7>l + (1 + l)ul] u z (3) 

which is a system of ordinary differential equations fully 
characterizing the transverse profile of the considered 
nonlinear guided waves. Without loss of generality we 
consider solutions of Eqs.([3]) with definite parity where 
u x and u z are spatially even (u x (^) = u x {— £)) and 
odd (w.z(£) = — u z (— £)), respectively and, as a con- 
sequence, we adopt the boundary conditions u x (0) — 
u x o, u z (0) = and u x (+oo) = u xoo , u z (+oo) = 
u zoo . Since u x (^) and it z (£) have to asymptotically 
approach two constant values, their first and second 
derivative vanish for £ — > +00 so that, exploiting the 
boundary conditions, we require the right hand sides 




FIG. 2: (a) Profiles of the z- component of the Poynting vector 

(see Eq.©) normalized with So = <J (eoe 3 )/(4/ioMX 2 ) corre- 
sponding to the fields reported in Fig.l. Each profile is char- 
acterized by an off-center positive part (black portion) and a 
central negative part (red portion), (b) Plot of the field S/So 
(arrows) in the plane (£, Q) corresponding to the nonlinear 
guided wave with u xoo = 0.65 of Fig.(f). The color is related 
to the local value of S z /So. Note the reversing of S along the 
transverse £ axis. 



Therefore we obtain j3 — y/2~/(l — 2u xoo )/{\ + 7) and 



u zoo = y[I — (1 — j)u xco \/(l + 7) from which we note 
that u xoo < 1/2 is a necessary condition for the existence 
of the considered nonlinear waves. In order to prove their 
existence, we exploit the fact that the system of Eqs.Q 
is integrable [l5[ since it admits the first integral 

F{u x ,u z ) = {(3 2 - l)u 2 x - u 2 z + _(«*+«*) + (! - 7 ) U 2 U | + 



■1 [(/3 2 - !) + (! + j)u 2 x + (1-J)u: 



l 2 2 



(4) 



or, in other words, the relation ^F{u x (^), u z (£)) = 
holds for any solution u x (£), u z(0 °f Eqs.Q. Evidently, 
after substituting the obtained (3 and u zoc into Eq.((4]), 
F has a stationary point at (u xoo ,u zoo ) and the guided 
waves are represented by curves of constat F in the plane 
(u x ,u z ) joining (u x q,0) to the stationary point. There- 
fore, requiring that F has a saddle point at {u xoo ,u ZOQ ) 
and exploiting the above necessary condition, we con- 
clude that the considered nonlinear waves exist in the 
range 



1 



(7+1) 5 
7 



< u x 



< \ - 



(5) 



of Eqs.Q to vanish at 



and 



which is always not empty since 7 > 0. In addition the 
relation F(u x q, 0) = F(u xoo , u zoo ) yields the possible val- 
ues of u x (0) — u x q corresponding to the asymptotical 
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the form E(x, z) 



[E x (x, z)e x + iE z (x, z)e z ] since, if 



FIG. 3: Metamaterial layered structure able to support trans- 
verse power flow reversing of TM fields, consisting of alter- 
nating slabs of a negative permittivity dielectric (ND) and a 
standard nonlinear cubic medium (NL). 



value u x (+oo) = u xoo . In Fig.l we plot the profiles of u x 
and u z corresponding to different values of u XOQl spanning 
the range of Eq.©, for 7 = 0.5 obtained by numerically 
integrating Eqs.Q with the above boundary conditions. 
The power flow carried by these waves is described by the 
time-average Poynting vector S = ^Re(E x H*) which, 
exploiting Eqs.© and the first of Eq.©, becomes 



e e° 



1 



-{1- [(ul + u 2 z )+ 7 (ul-u 2 z )]}ult 



(6) 

i.e., for the considered waves, is purely along the z— axis. 
In Fig. 2 (a) we plot the profiles of S z evaluated for the 
fields reported in Fig.l, and we note that the sign of 
S z is not constant along the transverse profiles, a region 
where S z < (red portion of the curves) existing around 
£ = 0. This reversing of the power flow along the trans- 
verse profile of the nonlinear guided waves is particularly 
evident from Fig. 2(b) where we draw the vector field S 
on the plane (£,C) for one of the fields of Fig.f. In or- 
der to physically grasp and discuss this unusual effect we 
recast Eq.© in the form 



S = 



2/3 ^ 



-.d x e: 



(J) 



where use of Eqs.([2]) and fTJ has been made, from which 
it is evident that the transverse power flow reversing is 
a consequence of the sign flipping of D x along the wave 
profile while E x does not change its sign. This implies 
that, regardless the absolute sign of the fields, the over- 
all effective dielectric response undergoes a sign reversing 
due to the fact that, in the first of Eq.([T]), the nonlin- 
ear cubic term can be both smaller and greater than the 
linear part, depending on the local field strengths. It 
is worth noting that, although we have discussed this 
effect using the considered nonlinear guided waves ad- 
mitting analytical treatment, the phenomenon is more 
general and holds for the wider class of TM fields of 



\d x E x \ < k\E x \ and \d x E z \ <C k\E z \ (i.e. the field manly 
propagates along the z— axis) it is simple to obtain from 
Maxwell equations that S z = ^Re[D*(x, z)E x (x, z)] 
and the transverse power flow reversing can take place 
through the just discussed mechanism. We conclude that 
the predicted power flow reversing is a signature of the 
extreme nonlinear regime where the cubic nonlinear con- 
tribution to the medium polarizability is not a mere per- 
turbation of the linear part. In this sense the medium 
behaves as a metamaterial whose character (positive or 
negative dielectric constant) locally depends on the field 
intensity. It is worth stressing that the discussed power 
flow reversing is very different from the effect that, in left 
handed metamaterials, the Poynting vector is antiparallcl 
to the carrier wave vector which is a consequence of the 
fact that, in such media, e < and /1 < (with n < 0). 
On the other hand, in our case, /j, > and the sign of the 
power flow is not uniform being controlled through the 
field intensity. Note that such an extreme condition can 
be achieved when the field intensity \E\ 2 is comparable 
or greater than e/x, so that, in standard materials where 
e is generally of the order of unity and x is very small (of 



the order of f m /V in semiconductors 



ery sr. 



the re- 



quired intensity is so large to rule out the whole discussed 
phenomenology. However, if a metamaterial is employed 
where e can be chosen to be much smaller than unity, 
the intensity threshold can be reduced to the point of 
making the extreme nonlinear regime accessible even for 
intensities much smaller than those employed in standard 
nonlinear optics experiments. 

Even though the use of a metamaterial (e <C 1) 
makes feasible intensities able to trigger the above linear- 
nonlinear competition, the main issue remains of finding 
a medium whose dielectric response is, in the considered 
intensity range, purely cubic. In fact, the first of Eqs.fT]) 
is a power series expansion of the constitutive relation 
D = D(E) in the field strength E and therefore, if the 
third order is comparable with the first one, one generally 
has to consider higher order terms. In order to show that 
the discussed extreme nonlinear regime can effectively be 
achieved, consider the metamaterial structure reported in 
Fig. 3 consisting of alternate linear metamaterial and non- 
linear medium layers of thickness d\ and d 2 respectively. 
The metamaterial is a negative dielectric (ND) whose 
constitutive relation is D = eoeiE (ei < 0) whereas the 
nonlinear medium (NL) is characterized by the constitu- 
tive relation D = e e 2 E - e X2[(E • E*)E + 7 (E ■ E)E*], 
i.e. it is a standard isotropic defocusing (e 2 > 0,X2 > 0) 
Kerr medium. The media relative permeability are /ii 
and /i2, respectively. If the spatial period d\ + d 2 is much 
smaller than the field vacuum wavelength 2itc/uj, the TM 
field propagating through the structure experiences the 
effective response described by Eqs.([T]) and characterized 
by the spatially averaged parameters 

e = /ei + (l-/)e 2j 
X = (1 - f)X2, 
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where / = d\/(di + cfe) is the fraction of negative di- 
electric. From these relations, it is evident that suitable 
values of ei and / can be chosen so that < e <C 1, i.e. 
the overall medium effective response coincides with the 
one considered in present Letter. Most importantly, the 
medium is able to support the extreme nonlinear regime 
since if the field is such that \E\ 2 ~ e/x, at the same 
time one has that \E\ 2 <C ti/x2- Therefore the nonlin- 
ear medium layers (NL) are in the presence of a field 
for which their response is purely cubic and, as a conse- 
quence, the overall averaged structure response is purely 
cubic as well. 

Even though the physical origin of the extreme non 
linear regime is entirely dielectric, magnetic properties of 
the metamaterial can play an important role. In fact, 
from Eqs.([2]), Eq.© and the relations £ = ^/efi(uJ / 'c)x , 
C = ^/ejl(uj / c)z, we note that fi plays the role of a scal- 
ing parameter for both the field physical sizes and inten- 
sity. This implies that, the field propagating through a 
medium with fi > 1 can be obtained from the field prop- 
agating through a medium with fj, = 1 by shrinking its 
sizes and reducing its intensity by a factor l/-»//L As a 
consequence the intensity required for observing the ex- 
treme nonlinear regime can be further reduced exploiting 
the range of available permeability fj,, achievable through 
specially designed metamaterials. 

As a specific example, consider a TM field of wave- 
length A = 1.064/im propagating through a layered meta- 
material structure for which e\ = —51.991 (coinciding, 
for example, with the real part of the silver permittivity 
at the considered A [13), /ii = 1 and £2 = 12.041, fi2 = 1, 
X2 = 8.742 x 10- 19 m 2 /V 2 and 7 = 0.5 (which are the 
linear and nonlinear parameters characterizing the GaAs 
at the considered wavelength 

[Mill)- From Eqs.© we 
obtain, for / — 0.188, the average effective permittivity 
e = 0.003 and the average permeability fi = 1. Consider 
now the nonlinear guided wave whose power flow is re- 
ported in Fig.(2)b which is characterized by u xoc = 0.65 



and a transverse dimensionless width A£ ~ 2 (A£ also 
coincides with the width of the transverse portion of the 
field where the Poynting vector is antiparallel to the prop- 
agation direction). The physical width of the considered 
wave is Ax — (A/27r)(A£/-y/e) ~ 6.2/im whereas the max- 
imum of its normalized Poynting vector is S z /Sq ~ 0.1 
(see Fig. 2) so that the wave is characterized by the inten- 
sity S z = 0.1 So ~ 3 MW/cm? . It is worth stressing that 
the considered micron-sized confined wave is observable 
with an intensity (~ MW/cm 2 ) much smaller than that 
(~ GW/cm 2 ) required for exciting a spatial soliton (of 
the same width and at the same wavelength) propagat- 
ing through a GaAs sample Referring to the above 
example, from a conceptual point of view, if everything 
is left unchanged but we choose the extreme permeabil- 
ity fi\ = —0.234, the considered nonlinear guided wave 
is characterized by the (sub-wavelength) physical width 
Aa; ~ 0.57 [im and the intensity S z ~ 0.27 MW/cm 2 , i.e. 
they are both reduced roughly by a factor of 10. 

In conclusion, we have shown that nonlinear meta- 
material with very small dielectric permittivity can sup- 
port the propagation of electromagnetic beams exhibiting 
transverse power flow reversing, i.e. the Poynting vec- 
tor changes sign along their transverse profile. Such an 
unusual phenomenon is one of the manifestations of the 
underlying extreme nonlinear regime where the nonlinear 
contribution to the polarizability can even exceed the lin- 
ear contribution, a situation never occurring in general 
nonlinear optical setups. It is worth stressing that the 
novel extreme nonlinear regime and its foreseen manifes- 
tations are observable even at relatively small intensities 
since the metamaterial has a small linear dielectric per- 
mittivity. In addition, we have noted that metamaterials 
with relevant magnetic properties can allow the observa- 
tion of the novel regime even for very small intensities 
and characterized by a very tight sub-wavelength field 
confinement. Therefore we conclude that the designing 
of complex metamaterials where dielectric nonlinearity 
combines with a pronounced magnetic response, can play 
a fundamental role for conceiving sub- wavelength nonlin- 
ear devices. 
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